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SUMMARY 


The present understanding of shock-layer radiation in the low density regime, as appropri- 
ate to hypersonic vehicles, is surveyed. Based on the relative importance of electron excitation 
and radiation transport, the hypersonic flows are divided into three groups: weakly ionized, 
moderately ionized, and highly ionized flows. In the light of this division, the existing labo- 
ratory and flight data arc scrutinized. Finally, an assessment of the nonequilibrium radiation 
computation methods for the three regimes in hypersonic flows is presented. The assessment 
is conducted by comparing experimental data against the values predicted by the physical 
model. 


NOMENCLATURE 


Aij Einstein transition probability for spontaneous emission 

B\ Planck function 

c velocity of light, cm/sec 

< c > average thermal speed, cm/sec 

C constant, equation 2 

e v average vibrational energy per unit mass, J/kg 

e v E average vibrational energy per unit mass at equilibrium, J/kg 

en average rotational energy per unit mass, J/kg 

cre average rotational energy per unit mass at equilibrium, J/kg 

h Planck constant 

k Boltzmann constant 

K rate constant, cm 3 /sec 

k absorption coefficient, cm -1 

n constant, equation 2 

N number density, cm -3 

p pressure, atmosphere 

T translational temperature, K 

T a geometrical mean temperature, K 

Td characteristic temperature of dissociation, K 

T e electron temperature, K 

T ex excitation temperature (eq. 1), K 

Tr rotational temperature, K 

T v vibrational temperature, K 

v vibrational quantum number 

tl relaxation time deduced from the Landau-Teller equation, sec 



I. INTRODUCTION 


The physics of the compressed gas layer behind a shock wave surrounding a hypersonic 
object flying through a planetary atmosphere have been studied intensely since the early 
1960’s (refs. 1-18)., The intense heat generated as a result of the compression is distributed 
in the shock layer volume through a thcrmochemical relaxation process. The translational 
temperature immediately behind the shock may be as high as 20,000 to 60,000 K. Due to the 
high velocity of the vehicle the characteristic flow time through the shock layer becomes com- 
parable to the relaxation time constants of these thermochemical processes and, as a result, a 
major volume of the shock layer is in a state of thermochemical nonequilibrium. The radiation 
from the volume in thcrmochemical nonequilibrium is enhanced, about 2-15 times higher than 
its equilibrium value (ref. 16). In such an environment, the nonequilibrium radiation plays a 
larger role in the radiative heating of the vehicle body. In order to properly design the heat 
shields of such vehicles one must be able to compute the radiative and convective heat load. 
The convective heat transfer rates can be estimated with confidence provided the surface cat- 
alytic reactivity is known (ref. 19). The nonequilibrium radiation phenomena is still not well 
understood, especially, at high altitudes (> 50 km) and at velocities above 10.0 km/sec and 
below 6.0 km/scc. 

There are three basic features which arc unique to nonequilibrium radiation: 


1. For a given transition, the absorption coefficient, «, and the emission coefficient, c must 
satisfy the principle of detailed balance unique to that particular transition, independently 
of all other transitions occurring concurrently 


i = KB\{T ejc ), B\(T ex ) 


8A he exp(—hc/XkT ex ) 
A 3 * 5 1 — exp(—hc/XkT ex ) 


( 1 ) 


where B\(T CX ) is the Planck function dictated by the excitation temperature for a par- 
ticular transition. This condition is satisfied, whether the gas is in radiative equilibrium 
or not. If radiative equilibrium exits, the total absorption coefficient, K t =Y1 K ar *d the 
total emission coefficient, (t ~ arc also related by the same excitation temperature. 
In nonequilibrium radiation a single excitation temperature docs not exist. 


2. During the relaxation process a state of chemical nonequilibrium exits. The populations of 
individual radiating species arc functions of time (or distance), and must be determined by 
solving a set of chemical reaction rate equations simultaneously with the thermodynamic 
and gasdynamic equations. 

3. The internal energy distribution for the molecules as well as the atoms can not be uniquely 

defined. The molecular populations in individual rotational states is defined by a ro- 

tational temperature, Tr. The distribution electron energy is defined by an electron 
temperature, T r . The populations in atomic and molecular electronic levels can not be 
defined by a single temperature. Instead, they arc determined by solving another set of 
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rate equations under the so-called quasi-steady-state (QSS) approximation (ref. 7). Re- 
cent computations (ref. 20) suggest that the vibrational populations of the molecules can 
not be determined by assuming a single vibrational temperature, T v and, in principle at 
least, should be determined by solving a set of master equations (ref. 20). Even in the 
simplest models, one must define a vibrational temperature for the vibrational degree 
of freedom of the molecules. Also individual molecular species may not have the same 
vibrational temperature (ref. 21). 

During the Apollo era to understand the basic physics of nonequilibrium radiation and 
compute the radiative load to the vehicle, the radiative phenomena were studied extensively 
through laboratory (refs. 1-3) and flight tests (refs. 10-14). In a concerted effort in the labo- 
ratory, shock tube and ballistic range tests were conducted along with theoretical calculations. 
Shock tubes were operated (refs. 1-3) at shock velocities up to 10 km/sec. There was a lower 
limit on the gas density in these tests because of intrinsic limitations of a shock tube, such as 
boundary layer growth and short test times. At the shock velocity of 7 km/sec, the lowest ini- 
tial pressure ahead of the shock wave was 0.05 Torr (~70 km altitude). At the shock velocity 
of 10 km/scc, the lowest pressure was 0.1 Torr (~65 km altitude). One of the major findings 
of the laboratory data was that the luminosity of air is zero immediately behind the shock, 
rises to a peak, then, decays to a smaller value and achieves a plateau (fig. 1). This plateau 
suggests the onset of a steady state (very close to equilibrium). This nonequilibrium radiation 
overshoot was interpreted as being due to the thermochemical nonequilibrium. The molecules 
with high binding energy (~9.8 cV) require a long time to dissociate, as the molecule cascaded 
towards the dissociation limit by way of collisional energy transfers. Before it could dissociate 
to its equilibrium value, the temperature rose to a higher level than its equilibrium value, 
causing the luminosity overshoot. However, the first attempts to theoretically simulate the 
radiation overshoot were unsuccessful. The model should allow the reaction rates to be small 
at the initial stage of the relaxation process and, then, to grow exponentially near the final 
stages of the process. If one uses the conventional formulation for reaction rate, K, 

Reaction rate, K = CT n exp(—Td/T) (2) 

with the translational temperature as the controlling temperature, then the maximum radia- 
tion occurs immediately behind the shock, where the temperature is maximum. Attempts to 
attribute the incubation period to the nonequilibrium vibrational excitation have also failed, 
because peak in radiation occurs much later in time than predicted by the vibrational relax- 
ation theory using Landau- Teller formulation (ref. 22). Over the years several investigators 
(refs. 23-28) have suggested that the dissociation and vibrational relaxation processes be 
coupled together in order to explain this long incubation period to the peak radiation. 

Recently, Park (ref. 29), was able to theoretically simulate the slow relaxation and radi- 
ation overshoot for shock velocities in the range of 6-10 km/sec using a multi - tem perature 
model (T v = T e ^T R = T), by suggesting that an average temperature, T a = yJTT v controls 
the reaction rates in equation (2). The model also suggests that at very high temperatures 
the vibrational relaxation changes its character from Landau-Teller type to a slower diffusion 
type (ref. 29). 
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Since chemical and thermodynamic processes in the relaxation zone are binary, i.e. they 
result from two body collisions, the intensity of radiation from such reactions is proportional 
to the density and the thickness of the relaxation zone is inversely proportional to the den- 
sity. Consequently, as long as the binary mechanism dominates, the integrated radiation flux 
emitted from such a nonequilibrium region remains independent of density at a given velocity 
(refs. 2 and 3). This property is called the “binary scaling” law of nonequilibrium radiation. 
There are four mechanisms which may depress, either the nonequilibrium or equilibrium ra- 
diative load to a body, and thus violate the binary scaling law: 1) Truncation occurs (ref. 2) 
when the width of the nonequilibrium region is thicker than the shock detachment distance, 
2) collisional limiting (ref. 8) reduces the radiation intensity at low densities when the colli- 
sions to maintain the population of excited states against the depletion by emission are not 
sufficient, 3) radiation cooling occurs when energy loss by radiation is significant, and 4) the 
gas is not always optically thin and reabsorption in the gas may occur. 

Electrons play a dominant role in the excitation and emission process. Electrons are 
about a million times more efficient than neutral particles and few hundred times than ions in 
causing electronic excitations of atoms and molecules by collisions. Furthermore, experiments 
show (ref. 30) that high electronic states, which primarily contribute to the emission, remain 
in equilibrium with the electrons described by a local Saha equation. Thus, for singly ionized 
species ( N e ~ Ni ons ), the population of excited electronic states is maintained proportional 
to the square of the electron number density. In this case 

Spontaneous radiative power = hvAi i jN l g l ~ huAi^N^ (3) 

where Aij is the Einstein transition probability for spontaneous emission for the i —> j tran- 
sition, Ni is the number density of the ith electronic level and gi is the degeneracy of the ith 
level. The electrons not only affect the line radiation but also the continuum radiation as 
well. The continuum radiation, either due to free- free transition (oc N%) or due to bound- free 
recombination process (oc N e ) or due to density packed bound-bound line radiations (oc JVg), 
can add up to a significant portion of the total radiative power, especially in an optically thick 
environment. 

Thus it is seen that the nonequilibrium radiation is affected not only by the thermo- 
chemical nonequilibrium in the media, but also by the vehicle size and flight altitude. Also, to 
a certain degree, the nature of radiation depends upon the electron density. Todate, various 
radiation models have been reported in the literature, each of which capture some of the 
important features of nonequilibrium radiation. This survey paper examines the capabilities 
of these models and summarizes the status of our present day understanding of the physics 
of nonequilibrium radiation. As an example, the results obtained from the computer program 
named NonEquilibrium Air Radiation (NEQAIR) written by Park (ref. 7) at NASA Ames 
Research Center, Moffett Field, California, U.S.A., will be referred most frequently. The 
available experimental data, from laboratories as well from flight experiments is surveyed. 
Finally some of the results obtained using these models are discussed and compared with the 
experimental data. The nonequilibrium radiation as a whole is divided into three groups: 
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(i) Radiation at low ionization, N e /N < 0.01%. 

(ii) Radiation at moderate ionization, 3% > N e /N > 0.01%. 

(iii) Radiation at high degree of ionization, N c /N > 3%. 


II. PHYSICS OF NONEQUILIBRIUM RADIATION 


A. Atomic Radiation 


In order to calculate the total radiation from a bound-bound optical transition ( i —* j) 
one must compute the spontaneous radiation (oc JVj), the stimulated radiation (oc Ni), and 
the absorbed radiation (oc Nj). This requires the knowledge of number densities, Ni and Nj. 
As mentioned before N{ and Nj are related by a local excitation temperature (eq. (1)). Under 
nonequilibrium conditions the number density of any given state, ‘i’ of an atom or molecule 
changes with time and is a sum of the diffusion rate and the net chemical production rate. 
The chemical processes consist of collisional and radiative transitions. For internal states, 
however, the chemical production as well as chemical removal rates are usually much larger 
than the diffusion term. This leads to the QSS approximation 
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.. . . 
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( 4 ) 


here N e is the electron number density (the colliding particle). K(i,j) is the rate coefficient 
for i — > j electronic transition, and K(i,c) and K(c,i ) are the rate coefficients for i — > free 
and free — > i transitions, respectively. The rate coefficients are obtained by the integration 
of corresponding cross sections over collision energy from the threshold value to infinity. In 
general the colliding particles can be a heavy particle and an electron. However, since the 
electrons are several orders of magnitudes more effective in causing electronic or vibrational 
transitions than heavy particles, for these transitions the contributions of heavy particles may 
be neglected (ref. 32). Determining the number densities, Ni, by solving the QSS Equations 
(cq. (4)), in essence, we assume a multi electronic temperature model for each species while 
computing the radiative emission (and/or absorption). 
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In general, of the three atomic radiative processes, namely: 1) line radiation, 2) free-bound 
radiation, and 3) free-free radiation, the radiative intensity from the first two are affected by 
the populations of all internal electronic states. The free-free radiation is not affected by the 
internal state populations. For the computation of the line atomic radiation and free-bound 
continuum the populations of all electronic states are determined by solving a set of equations 
like equation (4). The transitions that need to be considered in eq. (4), can be divided into 
four groups: 

1. Transitions among high states (i > 4). A QSS condition is assumed to exist among 
these states. The rate coefficients for the radiative transitions can be calculated by 
using references 33 and 34. The collisional rate coefficients are derived by Bates 
(ref. 35) using the semi-classical theory of Gryzinsky (refs. 36 and 37) for a hydrogen 
force field. The collisional rates of Bates can be used for N and O with an effective 
quantum number i.c. replacing En/n 2 by (E 0 Q — E(i)). Here Eh is the ionization 
potential for hydrogen and n its principle quantum number. 

2. Transition among low lying states (i < 3). A QSS state does not exist among these 
states. However, they maintain a state of equilibrium among themselves and their 
populations can be described by a Boltzmann distribution. In order to compensate 
for the large departure from a hydrogenic structure for the three low lying states 
in atoms such as nitrogen, a correction factor can be used (sec ref. 38). Within the 
framework of Gryzinski’s theory this correction factor applies only to the initial states 
and is independent of the final state. However, at low densities and highly nonequilib- 
rium regimes the correction factor for a singlet-to-singlet transitions may be different 
than for a singlct-to-multilet transition. Thus the applicability of Gryzinski’s theory 
to the low density high velocity regimes needs to be further researched. 

3. Transitions from the ground state to the high excited state. Experimental data for 
collisional excitation are available (refs. 38 40). 

4. Transitions from the low lying states ( i < 3) to the high states. Due to lack of 
additional data the cross section for this transition can be assumed to be the same 
as for the transition from ground to the high excited states. 

In a collision dominated plasma, when the contribution of radiative transition to equation 
(4) is not significant, the population distribution of the electronic states becomes independent 
of electron number density and insensitive to the excitation rate coefficients, or to errors in 
their estimates. A collision-dominant condition exits in one of the following situations: 


a. Optically thick gas. 

b. Optically thin gas with sufficiently high number of electrons. 

A 

) » —M > j and i,j > 4 


( 5 ) 


6 


For nitrogen this conditions is satisfied for N e > 10 18 cm 3 . 

c. Partially optically thick gas. For a given optical depth there exists an upper limit 
for the electron number density (< 10 18 cm -3 ) when a partially thick gas becomes 
collisionally dominated. For example at T = 12000 K in nitrogen this limit is N e > 
10 16 cm' 3 . If the far ultraviolet radiation is cut off, then the limit for nitrogen is 
N e > 10 14 cra~ 3 . 

Of course at electron densities below the limiting value, the radiative transitions must 
also be included in the calculations. 

The integrated intensity from spontaneous emission is distributed over a spectral line 
shape. Most realistic line shape can be approximated by the Voigt line shape (ref. 41), 
which combines the Lorentzian and Gaussian line shapes into one. The Gaussian shape is 
produced by the thermal motion of atoms. Overall line shape, however, is due to combined 
effects of natural broadening, Stark broadening (electron collisions), Van der Waals broadening 
(other species collisions) and resonance broadening (same species collisions). At low densities 
of present interest, the line shapes can be assumed to be symmetrical about the centerline 
(ref. 41) and unshifted by collisions. However, in certain cases, asymmetric line profiles due 
to Van der Waals broadening may occur (ref. 42). 

When the electron number density is in excess of 3% (>9.0 km/scc) Kramers’ radiation, 
namely: free-free, produced by deceleration of an electron in the electric field of an ion or 
atom, and free-bound, produced by capture of an electron by an ion, also becomes important. 


B. Molecular Radiation 

The molecular spectra for any given electronic transition consists of a vibrational band 
system with a rotational structure. In order to compute the radiation from each of the rota- 
tional lines one must determine the rotational populations in each of the involved vibrational 
and electronic states. The excitation energies of molecular rotation are usually very small. For 
example, the temperature equivalent for oxygen is 2.1 K and for nitrogen is 2.9 K. Therefore, 
at temperatures of our interest, the quantum effects of molecular rotation arc not pronounced. 
The rotational levels tend to establish an equilibrium among themselves very quickly, for ex- 
ample, within 20 collisions in nitrogen, and their distribution can be described by a rotational 
temperature, T/j. Because of a “near-classical” nature of the rotation, a vigorous exchange 
between the translational and rotational degree of freedom occur during collisions. As a re- 
sult, the rotational temperature tends to follow the translational temperature. However there 
arc experimental evidences (refs. 4 and 5) indicating that the rotational temperature may lag 
behind the translational temperature for quite some distance in a flow behind a strong shock. 
Thus, for example in general one needs to define T/j and T separately. 

To first order, the nonequilibrium vibrational population may be described by a single 
vibrational temperature, T„. However, recent computations show that during relaxation the 
vibrational states, with the exccpl ion of lowest levels, maintain a high degree of noncquilibrium 
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among themselves and it is not possible to define a single vibrational temperature for all 
the vibrational levels. Therefore, whenever it is possible, one should use multi-vibrational 
temperature description of the gas medium. 

Once the vibrational temperature T v , the rotational temperature T«, the electron tem- 
perature, T e and the translational temperature, T are given, all that remains is to determine 
the populations in the electronic states. To do this, a set of equations similar to equation (4) 
arc solved. An effective rate coefficient K(i,j ) for i — » j electronic transition is defined by 
summing up the rate coefficient for one rotational-state-to-rotational-state transition, K(J, J') 
over all final vibrational and rotational states of ith electronic level and averaged over the ini- 
tial vibrational and rotational states of the jth electronic level (ref. 7). It is also assumed 
that the electrons are ineffective in causing rotational quantum changes and the cross section 
for the J — » J' transition is the product of the Frank-Condon factor and an electronic part, 
independent of either vibrational or rotational quantum numbers. The excitation cross sec- 
tion for air molecules has been measured (refs. 40 and 43-46) for T„ = T = T/* = 300 K for 
almost all transitions in which the initial state is the electronic ground state and are of the 
order of 10 _17 cm 2 , except for Nj , for which it is about 10~ 16 cm 2 . This room temperature 
data can be extended to any arbitrary combination of T, T„, Tjj and T e using the standard 
formulation (ref. 7). However, extending the room temperature data to the high velocity low 
density (U s > 10.0 km/sec, pi <0.1 Torr) by using these standard formulations may lead to 
errors. For transitions in which the ground state is not the initial state, the cross sections are 
not known and they arc assumed to be the same as the known cross sections with the same 
final state. 

For conditions where the electron mole fraction is at least 0.01%, the electrons play the 
dominant role in the excitation of electronic states (eq. (4)). However at the low end of electron 
densities, the excitation by heavy-particle impact may be significant. 

Some electronic states, such as the C 3 n u state of N 2 (upper state of N 2 ( 2 + ) band system) 
and the B 3 T,~ state of O 2 (upper state of Schumann-Runge Band System), need special treat- 
ment. They lie above the dissociation limit. Because they are so high, collisional excitation 
from the lower electronic states is unlikely. Instead, they are populated by collisions between 
atoms through a process called predissociation. Depopulation is by the reverse process of 
predissociation and by radiative decay. By applying the principle of detailed balance (ref. 7) 
one obtains the populations in the two states. As mentioned above, the populations of other 
states arc determined by solving the QSS equations in NEQAIR. 

From each electronic state of a molecule, dissociation by heavy particle impact occurs 
(electrons arc less efficient in this regard). Available rate data (ref. 47) can be utilized in 
deducing the rate coefficients for collisional dissociation from various electronic states. 

The radiative transition probabilities for molecules are also available in the literature 
(refs. 48 and 49). 
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III. NONEQUILIBRIUM CHEMISTRY 


A. Reaction Rates 

The gas medium, around hypersonic vehicles or in a hypersonic nozzle, most often, is a 
mixture of atoms and molecules. Under extreme temperature, they undergo a complex chain 
of chemical as well as thermodynamic interactions. As mentioned before, under nonequilib- 
rium the populations of individual species, which are numerous, remain functions of time (or 
distance). For synthetic air (78% N 2 , 22%0 2 ) computations, one has to consider about 11 
dominant species involved in about 22 different chemical and radiative processes and for a 
more realistic air composition (78% N 2 , 21% 0 2 , 0.93% A and 0.03% C0 2 ), one needs to 
consider 19 species with 38 processes. In order to compute the source term (net production 
rate) associated with the respective mass conservation equation for each species, one needs 
to know the reaction rates, as a rule, in a form similar to equation (2). The chemical rate 
data are available for only a few of these reactions from experiments; for the rest of the re- 
actions educated guesses arc made. Park (refs. 50 and 51), through an extensive study of 
flight and shock tube data, has compiled a set of chemical rate data (table 1, cq. (2)), which 
seems to simulate the thermochemical nonequilibrium in air, for shock velocities in the range 
of 6.0-10.0 km/sec (with moderate ionization), reasonably well. However, whenever the ther- 
modynamic conditions change above a certain level, such as, the degree of ionization increases 
or new species arc introduced (Martian atmosphere, for example) (ref. 52), the reaction rate 
data set invariably needs reexamination. 

Another point of major concern is the nature of controlling temperature in equation (2) 
to be used in high temperature nonequilibrium conditions. Park’s average temperature model 
(ref. 29) for dissociative reactions (table 1) seems to work well in predicting the nonequilibrium 
phenomena up to shock velocities of 10 km/sec. The model allows dissociation of a molecule 
by both centrifugal tearing, caused by rotational motion, and vibrational breaking. The 
square root of T (the model assumes Th = T) is a measure of the velocity associated with 
the rotational motion, whereas the square root of T„ is a measure of the velocity associated 
with the vibrational motion at the midpoint of its oscillations. The product conceptually 
represents an average energy level in phase space. In reactions where electrons are involved as 
the product of the reaction, it is recommended that T c is used as the controlling temperature, 
such as in ionization or recombination processes (table 1). 

In the absence of a sufficient number of electrons, the electronic excitation occurs pri- 
marily by collisions with heavy particles. Based on shock tube experiments, Lotz (ref. 53) has 
postulated that the cross section for electronic excitation by neutral species is zero at threshold 
and varies approximately as oc a.Log c (y)/y , where a is a constant on the order of 10 16 cm 2 and 
y is the ratio of the translational kinetic energy of the colliding particles to the threshold. For 
ionic molecules Lotz proposes that the cross-section is a constant (ref. 53). As an example of 
electronic excitation by heavy particles, the formation of excited NO molecules at low veloci- 
ties in air (refs. 47, 54, and 55) can be cited. The three body reactions N + 0 + M — > NO* + M 


9 


leads to the excitation of the A and B electronic states of NO, which are responsible for NO-7 
and NO-/? band systems. 


B. Multi- Temperature Model 

Wc have seen that while computing nonequilibrium radiation, one has to conduct a rigor- 
ous book keeping of all the electronic levels, possibly a score of vibrational levels, and locally 
defined rotational, vibrational and electron temperatures. For flow computations, such a rig- 
orous exercise for the internal energy quanta is prohibitively expensive and is not necessary. 
A simplified approach may be adequate. In this approach: 

1. In addition to translational energy, we recognize three separate pools of internal 
energy: 1) rotational energy, 2) vibrational energy, and 3) electronic energy. All the 
species contribute to these common pools of internal energy. In a more rigourous 
approach each species may be assigned to separate pools of energies, which in turn 
can be summed up to provide a global set of internal energies. 

2. Assign each pool of energy an associated distribution function and the respective 
temperature i.c. define the vibrational temperature T v , rotational temperature T/*, 
electron temperature T e , electronic temperature T e / ec t and translational tempera- 
ture T. 

3. Account for the interaction among the five pools of energy during the relaxation 
process. 

Rotational relaxation Due to the “classical” nature of the rotational energy, the 
rotational degree of freedom participates in a vigorous energy exchange with the translation 
mode. With that in mind, it is often customary to assume T = T/j. However, in a more 
accurate model one may define a Landau-Teller type relaxation for the rotational energy: 

degjjji) _ ejiEjT) - e^(Tfl) 
dt t lr 

where e R is rotational energy defined at a rotational temperature, T R . The rotational relax- 
ation parameter, t rr has been measured experimentally for molecules like hydrogen (refs. 56 
and 67). Data on air molecules are unavailable, primarily because the rotational tempera- 
ture follows the translational temperature very closely and experimentally it is difficult to 
distinguish between them. 

Vibrational relaxation The vibrational relaxation rate proposed by Millikan and 
White (ref. 58) shows that the Landau-Teller formulation is accurate at temperatures up to 
~5000 K, and holds reasonably good to approximately ~ 8000 K, although the equation tends 
to slightly underestimate the relaxation times at the higher temperature range. For higher 
temperatures the formula needs to be modified to take account of two phenomena: (ref. 7) 
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1) a limiting value of the cross section and 2) the diffusive nature of vibrational relaxation. 
The first modification arises from the fact that the Millikan and White correlation implies an 


unrealistically large cross section for vibrational relaxation at high temperatures. The correc- 


tion is introduced in the 
t' l is defined: 


form of a limiting cross section a v and an effective relaxation time, 


t' l =t l + t c , t c 


1 

< c > a v N 


where < c > is the thermal speed of the molecule and N the number density of the colliding 
particle. The numerical value of the a v is an unknown quantity and is adjusted to satisfy 
the existing experimental data. For air molecules, for example, a v = 10 -17 (50, 000/T) 2 gives 
satisfactory results (see fig. 2). 


The correction for the diffusive nature of the vibrational relaxation is suggested as follows: 
(ref. 59) 


de v (T v ) e vE (T) 


dt 


- e v (T v ) f Ti ~T V \ 

1 U-7W 


s-l 


(7) 


Landau— Teller correction factor 


where the subscript i refers to the initial condition immediately behind the shock wave. The 
parameter s is an increasing function of Tj, the maximum value being 3.5. A bridging formula 
is proposed in the following form 


s = 3.5 exp(— 5000/T s ) 

where T s is the translational-rotational temperature immediately behind the normal shock 
wave. 


A more accurate approach would be a bi-modal vibrational model. Theoretical inves- 
tigations (ref. 20) show that during nonequilibrium the vibrational population can not be 
described by a single temperature. The high vibrational levels tend to equilibrate among 
themselves and the free state very early during the relaxation, and are not in equilibrium with 
the other states. The lower levels follow the Landau-Teller type relaxation. This necessitates 
the introduction of a bi-modal vibrational model. This model may even account for the ex- 
cess energy, which still remains unaccounted for in the present day hypersonic flow equations 
(ref. 60). 

Electron and electronic temperature Electrons thermalize themselves very quickly; 
however, they do not exchange any appreciable energy by collision with heavy particles. There- 
fore, the electrons, whose distribution is described by an electron temperature, remain loosely 
coupled to the heavy species through vibrational and electronic excitation processes. As dis- 
cussed earlier, the electrons dominate the electronic excitation process, and if one assigns 
an electronic temperature to the excited electronic populations, T e j cct , then the assumption 
T e ~ T e i ec t may be justified. 

Early researchers suggested that a strong coupling exists between the electronic degree 
of freedom and the vibrational mode (ref. 61). This hypothesis is in accordance with the 
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fact that the radiation is quenched more effectively by molecular than by atomic collisions 
(ref. 32). Faizullov (ref. 62), while using the sodium line reversal technique, found that the 
populations of excited sodium atoms were strongly coupled with the vibrational temperature 
of the nitrogen molecule. There are other indications (ref. 63) that a strong coupling between 
the electronic and vibrational modes exists. 

The electrons arc also very efficient in causing vibrational energy transitions in molecules. 
In such electron-vibration (c-V) processes, multiple-level transitions also take place with signif- 
icant frequency, contrary to the Landau- Teller formulation. Lee (ref. 64), using the excitation 
and de-excitation rate coefficients based on computational quantum chemistry (ref. 65), theo- 
retically analyzed the role of electrons in the vibrational excitation of nitrogen molecules and 
has suggested the following expression for the time constant for the e-V process 

log(p e T e ) = 3.91 (logT e f - 30.3 6{logT e ) 

+ 48.90, 1000 K <T e < 7000 K 

log(p e r e ) = 1.30{logT e ) 2 - 9.09 (logT e ) ^ 

+ 5.58, 7000 K <T e < 50000# 

where p e is electron pressure and r e is the time constant. r e is to be used in the equation (7) 
in place of t' l for the c-V process. 

As it is clear from the above discussion that in order to simplify the computation, in 
certain cases, T e = T e i ect = T„ is acceptable. 


IV. COMPUTER CODES AND FLOW COMPUTATIONS 


Park’s (ref. 7) NEQAIR program incorporates most of the physics discussed in Section II 
using the basic structure of a code written by Arnold et al. (ref. 31) Two computer codes have 
been written to couple the one-dimensional fluid equations with the nonequilibrium radiation 
physics (NEQAIR): 1) STRAP (Shock Tube RAdiation Program) for the shock tube flow be- 
hind a normal shock, and 2) SPRAP (Stagnation Point RAdiation program) for the stagnation 
streamline in the shock layer over a blunt body. Park (refs. 59 and 66) has extensively used 
these codes to reproduce the experimental data (refs. 1-3, 9-15, and 67) taken in shock tubes, 
ballistic ranges, and flight experiments at shock velocities in the range of 6 10 km/sec. Park 
was able to reproduce most of the laboratory data, namely the 1) rotational temperature, 
2) electronic excitation temperature, 3) temporal variation of radiation intensities, 4) char- 
acteristic relaxation times, and 5) the ratio of nonequilibrium - to equilibrium-radiative heat 
fluxes (with the exception of vibrational temperature, which will be discussed later). In order 
to do that it was necessary to adjust the reaction rates, for which no experimental or analytical 
data arc available, by trial and error until a good agreement between the theory and experi- 
ment was achieved. A set of reaction rates for air constituents used in these computations is 
given in table 1. 
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Recently these program have been improved by Whiting (ref. 69) and combined into a 
single code. The improvements are: 1) generalization of chemical reaction schemes, 2) im- 
provement in the accuracy of partition functions and radiation calculation schemes, and 3) in- 
clusion of the species C, CO, CN, and C2 in the excitation calculations. All of these codes 
assume T eJec = T e =T V and = T in the fluid computations with the average temperature, 
T a = y/7% in the rate equations (eq. (2)). Several two- and three-dimensional fluid codes 
coupled with NEQAIR have been developed over the years. Two of the most notable of these 
are a program by Candler (ref. 70) and the program LAURA (Langley Aerothcrmodynamic 
Upwind Relaxation Algorithm) by Gnoffo (ref. 71). 

Chung Park (ref. 72) working on “Cassini” project extended the NEQAIR code to include 
methane, hydrogen and argon in order to simulate entry into Titans’s atmosphere. His fluid 
code used a three temperature model i.e. T e = T e i ec 7^ T v . Chung Park has also conducted 
emission radiation measurements in gas mixtures simulating Titan’s atmosphere and has com- 
pared them with his theoretical results. Nelson (ref. 73) also investigated the nonequilibrium 
radiation from Titans’s atmosphere by extending the SPRAP code to a host of hydrocarbons 
as constituents. 

For the simulation of low velocity rocket flights (~3.5 km/sec) Levin (ref. 54) has modified 
the NEQAIR code to include the electronic excitation of NO by heavy particle collisions and 
has compared the theoretical results with the flight data from the Bow Shock Ultraviolet 
Rockets (refs. 74 and 75) flights. 


V. COMPUTATIONAL RESULTS: DISCUSSION 


A. Air Radiation 

(i) Low ionization flows ( N e /N < 0.01%)- At low velocities (<5.0 km/sec), as dis- 
cussed in Section III, the electron population may not be great enough to dominate the elec- 
tronic and vibrational excitation process, and the heavy particle collisions play a major role. 
Since the excitation rates to discrete electronic levels of atoms and molecules by heavy par- 
ticle collisions are not very well known, modelling of nonequilibrium radiation in this regime 
is difficult. Strong coupling of vibrational energy with the electronic states of atoms and 
molecules (refs. 32 and 61 63) further complicate the physics of nonequilibrium radiation. In 
Levin et al’s computations (ref. 54) the electronic excitation of NO was modelled as a three 
body recombination process. The modified NEQAIR code was used to simulate the radiation 
recorded aboard the Bow Shock Ultraviolet Rocket flight experiment (ref. 74). The results 
show that the low altitude data («40 km, ~ equilibrium regime) agree with the computed 
results. However, the predicted nonequilibrium radiation using a multi-temperature model 
(T/t ^ T) was lower by as much as 4 orders of magnitude (at ~60 km). 

A possible explanation can be provided as follows. In recent shock tube experiments at 
Calspan (ref. 76), conducted at shock velocities in the 3.22-3.83 km/sec range, it was found 
that in the absence of NO, the vibrational temperature of N2 (infrared measurements at 5 /x) 
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shows a slow approach to its equilibrium value, in a manner NEQAIR predicts. In other 
words, the modelling of N 2 excitation in NEQAIR is reasonably accurate. As soon as NO 
forms in the flow, the vibrational temperature shoots up rapidly much faster than exhibited 
by N 2 vibrational temperature. It seems that NO molecules are probably produced in a 
vibrationally excited state. Possible mechanisms may be the formation of NO in reactions 
where the kinetic and chemical energies are directly converted into vibrational energy during 
the collisions. The vibrationally excited, ground electronic NO molecules may in turn interact 
to produce excited A-state molecules. This hypothesis seems to be in good agreement with the 
experimental findings about the strong couplings between vibrational and electronic states. 

This area of molecular excitation needs further work. Without a better understanding 
of these phenomena, our capability to compute nonequilibrium radiation at low velocities will 
not improve. 

(ii) Moderate ionization (3% > N e /N > 0.01%)- Most of the theoretical and exper- 
imental investigation on nonequilibrium radiation has been conducted in this regime. The 
low end, for example, may correspond to a shock velocity of 6.0 km/sec in 1 Torr air and the 
high end to a shock velocity of 10.0 km/scc into 0.1 Torr air. Spectral calculation of emitted 
nonequilibrium radiation is conducted by many investigators (refs. 1-5, 59 and 66) and com- 
pared with the available experimental data. The spectral data primarily have been recorded in 
shock tube experiments. The spectra at these conditions consists of NO-/3, N^l - ), N 2 (2 + ) 
and N 2 (l+) band systems. Typical equilibrium and nonequilibrium spectra recorded at a 
shock velocity of 10 km/sec in 0.1 Torr air (ref. 5) are shown in figures 3 and 4. As illustrated 
in figure 1 the nonequilibrium overshoot is very pronounced at these velocities. The ratio of 
peak nonequilibrium radiation to the equilibrium radiation (7 p // e , fig. 1 (c)), from various 
shock tube experiments arc plotted in figure 5. Four sets of the data: 1) at 4.75 km/sec in 
10 Torr N 2 , 2) at 5.54 in 3 Torr N 2 , 3) at 6.4 km/sec in 1 Torr N 2 , and 4) at 10 km/sec 
in 0.1 Torr air, were recorded by Allen et al 1 at AVCO Everett Research Laboratories, MA. 
The other two sets: 1) at 6.2 km/sec in 1 Torr N 2 , and at 10 km/sec in 0.1 Torr air, were 
recorded by the author and co-workcrs at NASA Ames (ref. 4 and 5). AVCO’s data were 
obtained by integrating the spectral emission data in the 5500-10000 A wavelength range and 
NASA Ames’s data were integrated in the 3100-7100 A range. The computed values for these 
conditions using STRAP’S two temperature model, T a = \/TT v for chemical nonequilibrium, 
also have been plotted in the figure. It should be remembered that for the computation of 
radiation STRAP uses NEQAIR, which models the electronic states of each species using a 
multi-temperature approach. The good agreement between the theory and the experimental 
data demonstrate that the physical model presented here simulates the radiation nonequilib- 
rium at moderate ionization with reasonable accuracy. It should be noted that the degree of 
nonequilibrium enhancement over its equilibrium value increases with increase in the shock 
velocity. At a certain shock velocity, most likely between 7 and 9 km/sec (pi ~ 0.1 Torr), 
Ip/1,, achieves a peak value and starts to drop. At 10 km/sec in 0.1 Torr air the ratio, I p /I e 
drops to R 2 4.5. It is important to note that at a given velocity the ratio, I p /I e is strongly 
dependent on the initial pressure before the shock. 
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As discussed before, the “binary scaling law” dictates that the nonequilibrium-to- 
equilibrium radiative heat flux defined as 

•k° 1 — (see fig. 1) 
h ■ to.l 

is a quantity independent of the initial pressure. Experimental values of this heat-flux ratio 
as function of shock velocity are compared with the computed values in figure 6. Once again 
a good agreement between the theory and the experiments is seen. 

The ability of the present model to compute detailed spectral radiation is put to test 
in figure 7, in which the experimental spectral emission data(refs. 1 and 5) at the point of 
peak intensity are compared with the synthetic spectrum produced by the STRAP code. The 
spectra are for a shock velocity of 10 km/sec in 0.1 Torr air. The synthetic spectrum is 
corrected for the instrument broadening by the input slit. The synthetic spectrum overesti- 
mates the emission in the UV (0.3-0. 5 //m) range and slightly underestimates in the infrared 
(0.55-1.5 /zm). The AVCO data show even lower emission in the UV range and higher emission 
in the infrared range when compared with the NASA Ames’s data. The overestimate of the 
UV emission, which is due to the Nj (1“) band system, could be caused by: 1) uncertainty 
of total Nj population estimates (» 50%), 2) effect of rotation on the dissociation process, 
and 3) the vibrational temperature may be too high. Recently, in temperature measurements 
at the point of peak radiation in nitrogen 4 and in air^ it was found that the rotational tem- 
perature was considerably lower than the translational temperature (fig. 8). This invalidates 
the assumptions T« = T made in STRAP and SPRAP codes and leads us to believe that a 
relaxation model for the rotational energy (eq. (6)) should be considered. 

The vibrational temperature of the species during the nonequilibrium is another matter 
of concern. In figure 9 the vibrational temperature measured by Allen (ref. 77) and by the 
Author 4 are plotted against the theoretical prediction for nitrogen (6.2 km/sec, 1 Torr). Two 
points should be made regarding this data. First, there is obviously a discrepancy between the 
AVCO and The NASA Ames’s data. Second, and perhaps more important, the vibrational 
temperatures measured at NASA Ames for N2(2‘ t_ ) and Nj(l ) are different in the same 
experiment. This result is not unexpected, since a single vibrational temperature does not 
necessarily exists for all species. Most likely the high vibrational states of N% ) tend to 

be in equilibrium with the free state (N% N + N + ). Park (ref. 50) determined 

a vibrational temperature, T* , describing the local equilibrium among high vibrational states 
and the free state, to be different from the ground state T„. T* computed this way follows 
the slow rise in the vibrational temperature measured by AVCO (fig. 9). To summarize, it 
can be stated that in general: 1) two different vibrational temperatures may exist, one for 
the ground states and another for the higher states, and 2) there may be different sets of two 
vibrational temperatures for each of the dominant species. 

The discussion of our present capabilities to compute radiation emission from high tem- 
perature air can be divided into two subgroups: 

a. Equilibrium radiation: The radiation plateau reached after the nonequilibrium over- 
shoot is most often believed to be caused by the onset of equilibrium in the gas. In fact, it 
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is due to the attainment of a steady state in which the populations in some of the higher 
electronic states may still remain much different from those of their equilibrium values. At 
low velocities (~6.0 km/sec) the radiative emission from this steady state is not very differ- 
ent from equilibrium air. However, at higher velocities and low densities (~10.0 km/sec) the 
difference becomes significant. It is interesting to note that this difference between the steady 
state and equilibrium emission is not uniformly distributed over the entire spectral range. Nj 
molecules contributing to the most prominent radiator, the First Negative Band System, seem 
to equilibrate near the plateau and hence the steady state and equilibrium emissions from this 
band system arc nearly the same. It should be noted here that CN violet radiation can signif- 
icantly alter the appearance of the (1~) band system, especially the (0,0) band system at 
3882.6 A and the (1,0) band system at 3582.1 A. Experiments (ref. 5) show that equilibrium 
radiation from 10 ppm CN impurity in the shock tube can be of the same magnitude as the 
emission from the Nj (l - ) band system. However, if the level of the impurities is properly 
accounted for, the agreement between the experiment and theory is excellent (ref. 78). 

On the other hand, at velocities higher than 10.0 km/sec and low densities, NEQAIR 
predicts much higher radiative emission from most of the atomic systems than measured 
experimentally. The discrepancy between the computations and the experiments does not seem 
be uniformly distributed over any particular atomic system, not even in the same multiplet 
system. For example, the computed emission from the 4p 2 * 4 5° — » 3 s 4 P nitrogen triplet system 
at 4151.46, 4143.42 and 4137.63 A falls, on average, within a factor of 5 of the experimental 
data. However, for the 4p 4 P° — * 3s 4 P multiplet system the computed values average about 
50 times higher than the experimental values. Individually, in the 4222.73, 4223.04 and 
4224.74 A triplet system, the computed values are about 61, 10.2 and 51 times higher than 
the experimental data, respectively (ref. 78). This large scatter in the discrepancies between 
theory and experiment persists even if a Boltzmann excitation is assumed. Two possible 
explanations are proposed: 

1. NEQAIR assumes that closely packed energy level achieve a local equilibrium during 
the relaxation process almost instantaneously and for that reason they are grouped 
together to form one pseudo state to simplify the QSS solution. It is possible that 
individual states of these multiplet systems, no matter how unlikely it may seem, 
populate and relax at different rates than their neighbors. If so, then one must 
formulate the QSS problem by accounting for all the relevant energy levels. 

2. The published ‘f’ numbers of the individual states used in NEQAIR may not be 
accurate. Most of the ‘f’ number data are obtained in a high current high pressure 
discharge where most spectral lines are broad. In a multiplet system, it is often not 

possible to distinguish between the adjacent spectral lines under these experimen- 
tal conditions. For simplicity, the contributions of individual lines in the multiplet 

system are assumed to be distributed according to the 1-s coupling. The emission 
experimental data (ref. 5) under low density high temperature conditions indicate 
that this may not be the case. Also ‘f’ number measurements are conducted under 
high electromagnetic forces. It is not clear how the ‘effective’ value of ‘f ’ numbers 
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will change, if any, under the relatively low ionization found in hypersonic flows. In 
any case, T numbers of the appropriate energy levels need to be investigated further. 

Also at velocities above 10 km/scc, the computed values of radiative emission from N 2 , 
are a few hundred times less than observed in the experiments. The under prediction from 
N 2 band systems makes sense in the following way. The molecular populations in certain 
electronic states may attain equilibrium with certain electronically excited free states (atoms). 
As a result, the code, on one hand overpredicts the atomic radiation and, at the same time, 
undcrpredicts the molecular radiation. This suggests that, in the experiments, the molecular 
populations in certain excited electronic states are lower than their equilibrium values and the 
atomic populations in certain the excited electronic states are higher. Although, at this stage 
the overall population of molecules and atoms are also in equilibrium at the temperature of 
equilibration (T = Tr = T v = T e ), experimental data suggest the distribution of electronic 
states still remains far from equilibrium. 

b. Nonequilibrium radiation: Emission in the nonequilibrium regime shows nearly the 
same behavior. At low velocities and higher densities (for example, at 6.0 km/scc in 1 Torr) the 
agreement between the experimental values and the prediction is reasonably good. At higher 
velocities and low densities (for example, at 10.2 km/sec in 0.1 Torr) the predicted emission 
from the band systems also agrees well with the experiments, provided the impurities 
like CN are properly accounted for. The level of impurities added in the computations to 
reconcile with the experimental data, however, are one to two order of magnitude lower than 
that required in the equilibrium case. 

As in the case of equilibrium radiation, at low densities and higher velocities the code 
under predicts the N 2 Second Positive and First Positive band systems and over predicts most 
of the N and O atomic systems. The code underpredicts the N 2 b§.nd radiation by about an 
order of magnitude and the overprediction of atomic radiation is about an order of magnitude 
less than in the equilibrium case. If all the collisions mechanisms are the same as in the 
equilibrium regime, radiative decay processes may be playing a dominant role in this regime 
and radiative life times may be comparable with the relaxation times. 

Beside the laboratory data, there are a host of flight data against which the nonequilibrium 
radiation model has been tested. Among the most notable flight data are: 1) Project FIRE 1 
(ref. 11) and FIRE 2 (ref. 12) data, 2) Planetary Experiment Test (PAET) (refs. 10 and 68) 
data, 3) Apollo 4 and 6 data (refs. 13 and 14). Flight experiments, unlike laboratory data, 
as a rule, provide data on bulk properties, such as integrated heat loads using calorimeters 
and radiometers. Apart from the flight altitude conditions (such as temperature, collision 
limiting etc.) the data is also affected by the phenomena caused by the vehicle’s geometry 
(such as truncation, ablation etc.). Also, since the vehicle flies in a real atmosphere, the trace 
elements found at the flight altitude may add additional complexities to the radiative heat 
transfer such as, observance of CN violet in PAET flights. Therefore, one must take these 
physical phenomena into account when analyzing the flight data. When one uses the flight 
data to validate or verify the physics of the nonequilibrium radiation, the added affects make 
the process more complex and the verification less convincing. At the same time, some of these 
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effects, such as collision limiting and truncation, provide another dimension in the verification 
process. 

The FIRE vehicles were approximately 1 /4-scale models of the Apollo command mod- 
ule and entered the atmosphere at velocities slightly higher than those of Apollo vehicles. 
The vehicles were made of beryllium and hence not ablating. The measured convective heat 
transfer agreed well with the available theories. The radiation data, at low altitudes, were 
in agreement with the equilibrium radiation computations. At higher altitudes, the radiative 
heating violated the “binary scaling law” and varied with diminishing density. However, if 
one takes “truncation” into account, it is possible to reproduce the flight data for a certain 
range of the flight (fig. 10). Due to the small size of the vehicle, many phenomena, such as 
viscous effects, finite shock thickness, small shock stand-off distance etc., play a role in the 
radiative transport. Therefore, only a portion of the flight data can be used to validate the 
physical model presented here. 

Over the years, attempts to reproduce the flight data mentioned here have been partially 
successful. The general consensus is that the true extent of nonequilibrium radiation in FIRE 
and Apollo experiments was not realized at the time. 

The preceding comparisons of experimental data and the theoretical prediction, in the 
moderately ionized regime, show that this model is fairly accurate in reproducing most of the 
experimental data. The model is reasonably accurate in reproducing integrated radiation in- 
tensities, and reproduces the spectral data at low velocities and high densities with acceptable 
accuracy. The model, however, does only a fair job of reproducing spectral data obtained in 
the ground tests at low density and high velocity. 

(iii) High degree of ionization (N e /N > 3%) The radiation emission in this regime 
is dominated by free-free transitions, bound-free transitions and strong atomic line radiation. 
Depending upon the optical path the atomic line tend to be self-absorbed in the media. Thus, 
in a moderately optically thick medium, which is a more realistic scenario for this regime, 
complex line by line spatial integration is required. In an optically thick approximation, of 
course, only free-free and bound-free radiation need to be considered. Therefore, emphasis of 
the model in this regime is more on the accuracy of radiative transport computations than on 
the chemical and/or thermodynamic nonequilibrium model. Because of high electron density, 
more charge transfer and associative ionization reactions arc added to the list of chemical 
reaction given in table 1 (ref. 52). 

In the highly ionized regime the nonequilibrium overshoots (end to decrease with increase 
in the shock velocity (fig. 5) leading to an undershoot. Such an undershoot can be seen in 
Wilson’s work (ref. 78), which is one of the very few experimental investigations available for 
this regime. With the selection of proper reaction rates (ref. 52), and by accounting for the 
alleged impurities in the experiments, the general trend of Wilson’s data can be reproduced 
qualitatively by the present model (fig. 11). Since the emission in this regime is proportional 
to N%, the plot in figure 11 is shown in terms of N%. 
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For the verification of computed spectral intensities, reliable and complete experimental 
data sets for this regime are not available. Depending upon the integration scheme in the 
radiative transport bookkeeping, discrepancies up to several orders of magnitude have been 
found in the spectral intensities. Even the integrated radiation heat flux, which under the 
“binary scaling law, ” should be density independent, has been found to scatter as much as 
one order of magnitude, depending on the method of computation. With a concerted effort in 
radiative transport integration techniques and the acquisition of reliable and complete spectral 
and integrated radiation data, our computational capabilities for this regime will be improved. 


B. Martian Atmosphere 

The Martian atmosphere consists of 95.7% CO 2 , 2.7% N 2 and 1.6% Ar (ref. 79). There- 
fore, a shock layer in the Martian atmosphere will contain C, CO, C 2 and CN, all of which 
have mechanisms for strong radiation. Of particular interest are the Fourth positive system of 
CO, the Red and Violet systems of CN, the Swan system of C 2 , and atomic lines of C and O. 
Reaction rates for various dissociation, ionization and vibrational relaxation of C0 2 , CO and 
N 2 are available in the literature (refs. 80-84). Fair amounts of theoretical and experimental 
data on the vibrational relaxation of various mixtures of C0 2 and CO have been made. The 
observed relaxation can be represented with reasonable accuracy by Millikan and White for- 
mula. A considerable amount of experimental data, valid for temperatures upto 15, 000 K, 
exist (refs. 82 and 83) on dissociation of C0 2 , CO and CN and on exchange reactions contain- 
ing these gases. Most of the high temperature experimental data are available for collisions 
with argon. The rates for collisions with other molecules and atoms are assumed (ref. 82) 
to be about 10 and 15 times higher than that measured in argon, respectively. Using these, 
tentative reaction rates applicable to Martian atmosphere can be compiled (sec ref. 52). 

The entry velocity into the Martian atmosphere range from 6-9 km/sec. No kinetic or 
radiation measurements have yet been made in the Martian gas mixtures. However, some 
theoretical and experimental investigations (ref. 85-91) have been made in various mixtures 
of C0 2 , CO and N 2 . These data can be used to assess the present model for calculation in 
the Martian atmosphere. In figure 12, the calculated equilibration distance (distance traveled 
in time to. 1 , fig- 1), reproduced here from reference 52, is shown and compared with the 
experimental data and other theoretical predictions for a post shock pressure of 0.15 atm. 
The prediction from reference 52 seems to agree with the experimental data of Arnold and 
Nicholls (ref. 86) approaching the mean value of Nealy’s data (ref. 88). Also it agrees with 
Howe’s (ref. 89) theoretical prediction. 

The computed values of radiative heat flux have been found to be very sensitive to line 
shapes of the transitions. In an optically thick environment, if one assumes a thin line, 
its center region, which contains most of the radiative energy, is self absorbed and the net 
radiative energy to the wall, contributed by the wings of the line, is bound to be smaller than 
in an optically thin gas. However, if one assumes a thick line, due to the smearing effect, the 
integrated area under the wings becomes appreciably larger than the thin line case. As a result, 
the computations estimate a larger heat flux to the wall under the thick line assumption. 
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Computations show that for flight velocity of 6-9 km/sec, the predicted equilibration 
distance varies between 9 and 27 cm. Thus, computations of the radiative heat load on an 
entry vehicle whose shock stand off distance is in the 9-27 cm range (or less), require the 
nonequilibrium radiation model. An equilibrium model will under predict the radiative heat 
load in those cases. 


C. Titan Atmosphere 

The atmosphere of Titan is known to be composed of molecular nitrogen, methane, hy- 
drogen and possibly argon. However, the specific mole fractions of the constituents are not 
well known (ref. 92). For the heat load analysis (ref. 72), a composition of 95.7% N2, 3.0% 
CH4, 1% argon and 0.3% hydrogen and an entry velocity in the range of 5-7 km/sec has been 
assumed. Under these conditions CN will readily be formed and is considered the dominant 
radiator along with minor contributions from C2 Swan band system. A list of chemical reac- 
tions and the reaction rates needed for the computations of radiative heating load in Titan’s 
atmosphere has been compiled by C.P. Park (ref. 72) and Nelson (ref. 73). C. P. Park also 
conducted emission measurements in a shock tube. He produced the CN emission data using 
a three temperature (T, T„ and T e ) and the average temperature T a = \JTT V model for the 
chemical processes and a modified NEQAIR code for the radiation calculations. 


CONCLUSIONS 


Based on our current understanding of the nonequilibrium phenomena, our computing 
capability, and the role of radiative transport and the thermochemistry in the physical model, 
we may divide the hypersonic flows into three groups: 1) low ionization flows, N e /N < 0.01%, 
2) moderate ionization flows, 3% < N e /N < 0.01%), and 3) highly ionized flows. 

Our basic understanding of nonequilibrium radiation in the low ionization flows is rather 
poor. The role of heavy particles in electronic excitation and the coupling between the vi- 
brational and electronic energy modes are, especially, not well understood. There is evidence 
that the published data on ‘f’ numbers for certain transitions may not be applicable to the 
regime of hypersonic flow. Also it seems that the QSS formulations performed for computing 
radiative emission (or absorption) must be extended to a larger number of energy levels. The 
experimental data being gathered at the E.A.S.T. Facility of NASA Ames Research Center 
have a wealth of information on these radiative processes. A thorough analysis of the data will 
improve our understanding of the relevant mcchnisms leading to the development of improved 
codes. 

The flows with moderate ionization are the most researched, and most understood. The 
spectral emission radiation for this regime can be computed with reasonable accuracy. The 
computed total integrated radiative heat flux from the shock layer for this regime agrees 
well with the shock tube, ballistic range and, in most cases, flight data. Other parameters, 
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such as the relaxation time and equilibration distance are also in good agreement with the 
experimental data. 

Our ability to compute the nonequilibrium radiation from highly ionized flows is also poor. 
The problem lies, not as much in the basic understanding of the physics involved, but in the 
complex bookkeeping of the radiative transport phenomena. Due to the significant degree 
of radiative cooling, radiative transport calculations must be done with improved accuracy. 
Also, reliable experimental data against which the computations can be verified, are limited 
in number. A low density shock tube facility is under development at the E.A.S.T. Facility of 
NASA Ames Research Center to address this problem. 

The computation of radiation from planetary atmospheres other than air is, in princi- 
ple, within our present day capabilities, since these flows fall into one of the last two cate- 
gories: moderately and highly ionized. Experimental data, however, are needed for verification 
purpose. 
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Table 1. Reaction rate coefficients for air; T r is the controlling temperature, C is in cm 3 mole -l s. -1 
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Figure 1. Concepts of nonequilibrium radiation: 1) shock layer, 2) approximate temperature 
profiles, and 3) nonequilibrium radiation characteristics. 
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Figure 2. Vibrational relaxation times as functions of temperature. 
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4500 5000 5500 6000 6500 

Wavelength, A 

Figure 3. Nonequilibrium radiation spectrum. Shock velocity = 10.20 km/sec, Pi = 0.1 Torr 
in nitrogen, image intensifier gate = 200 ns, spectrometer slit width = 60 /j.m, spectral reso- 
lution = 0.56 A/diode. Wavelength range 2800-4600 A. 




34 


in nitrogen, image intensifier gate = 1.0 fjsec, spectrometer slit width = 60 /jm, spectral 
resolution = 0.54 A/diode. Wavelength range 2800-4600 A. 



Figure 5. Ratio of peak nonequilibrium radiation to the equilibrium radiation as recorded in 
various shock tube experiments. 



Figure 6. Nonequilibrium-to-equilibrium radiation heat-flux ratio as function of shock velocity. 
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COMPARISON BETWEEN CALCULATED AND MEASURED 
SPECTRA AT PEAK-RADIATION POINT 

= 0.1 torr, V^j * 10 km/sec 



.2 .4 .6 .8 1-0 2-0 


Wavelength, nm 

Figure 7. Nonequilibrium spectral emission data compared with the synthetic spectrum. 



Distance (x, cm) 


Figure 8. The measured rotational temperature compared with the AVCO data (ref. 77) 
and with the theoretical values computed by Park (ref. 59) using a) average temperature, 
b) vibrational temperature, and c) translational temperature. 
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STRAP computations based on N 2 7 

STRAP computations based on Upper State 

O Allen 77 



Distance (x, cm) 


Figure 9. The measured vibrational temperature compared with the AVCO data (ref. 77) 
and with the theoretical values computed by Park (ref. 59) using a) average temperature, 
b) vibrational temperature, and c) translational temperature. 



Time from launch (sec) 

I I L I I _J 

90 80 70 60 50 40 
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Figure 10. Comparison between calculated and measured stagnation-point radiative heat 
fluxes for FIRE 2 (ref. 66). 
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• Wilson 78 
STRAP computation 82 



Distance from shock wave (cm) 


Figure 11. Comparison between Wilson’s experimental data and the data predicted by 
STRAP. 
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Figure 12. Comparison between the equilibrium distance at post shock pressure of 0.15 atm 
in Martian atmosphere predicted by the present model and the existing experimental data 
(refs. 85-88) and other theoretical predictions (refs. 89-91). 
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